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We find that the Measurement Based Quantum Computing (MBQC) search algorithm on an 
unsorted list is not the same as Grover's search algorithm (GSA). More importantly, the MBQC 
search algorithm is exponentially faster than both GSA and its classical counterpart. We develop 
and describe the numerical formalism we utilize that confirms the MBQC exponential speedup over 
the GSA and classical search algorithms. 
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PACS numbers: 03.67.Ac,03.67.Lx 

Introduction - L. Grover first showed that a speed up 
over the brute force classical search (CS) of an unsorted 
list was possible using quantum mechanics. His search al- 
gorithm, now known as Grover's search algorithm (GSA) 
[l| helped to spur interest in all manners of quantum com- 
puting due to the enhancement in computational speed 
over the corresponding classical algorithms. Grover's al- 
gorithm takes advantage of a process called amplitude 
amplification to act simultaneously on the full state vec- 
tor of quantum amplitudes []], Hf- It has been shown 
by C. Zalka [3] that the optimal speedup achievable by 
amplitude amplification is fl(V2™) where 2 n is the num- 
ber of elements in the unsorted list composed of n-bit 
strings (CS scales as ft(2 n ~ 1 )). However, the question 
remains if GSA is the optimal quantum search procedure 
for an unsorted list. It has been assumed that the Mea- 
surement Based Quantum Computing (MBQC) search 
method described and implemented on a photonic clus- 
ter state by Zeilinger et. al [4], is equivalent to GSA. 
We have found that this is not the case. It is therefore 
reasonable to expect that the MBQC search algorithm 
will have a different execution speed than GSA. We find 
that the MBQC search can be implemented exponentially 
faster than GSA and CS, with the MBQC search scaling 
as ft(log 2 2 n ) = Q(n). 

MBQC and its numerical formalism - MBQC operates 
by producing highly entangled multi-qubit states called 
cluster states. Computations are then preformed by a 
time ordered sequence of measurements on individual 
qubits in this entangled resource. We construct the clus- 
ter state by first forming the Kronecker product of N 
photons, each initialized to the state |+) = (|0)-h|l))/v / 2- 
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To complete the construction of the cluster state we 
create the entanglement resource between the qubits 
by preforming the CZ operation between nearest neigh- 
bors. The four CZ operations for the box cluster state 



(see Fig. QJi) are shown in (j2j), 
|*in>= II CZ iJ \i;) = CZ h2 CZ h3 CZ 2A CZ 3A \ 1 p). (2) 



Here we note that z, j is the label of the control and tar- 
get qubits respectively in the CZij operation which is a 
4x4 matrix with entries (— l) bi bj for the binary computa- 
tional basis entries bibj £ {0, 1}. In ((2]) nn is a list of all 
nearest neighbors. The CZij operations create the nxm 
rectangular cluster state resource the MBQC search re- 
quires with N = nm qubits. An algorithm, or any com- 



li Output 
4 



ll l 2 Output 

7 
8 



a) b) 
FIG. 1: Square cluster states. Two cluster states able 
to perform an unsorted search on a) 2 2 = 4 elements using 
\^box) and on b) 2 3 = 8 elements using |^C9)- Each circle is a 
single qubit initialized into the |+) state. The lines connecting 
neighboring qubits indicate the action of an entangling CZ 
gates. I\ and I2 are MBQC iterators on columns of qubits. 

putation, then proceeds as a series of single-qubit mea- 
surements in which the bases for the subsequent set of 
measurements is determined by the outcomes of the prior 
measurements. This can be considered an ordered series 
of Kronecker products of single qubit rotations yielding 
a factorizable unitary matrix. Equivalently, we may ro- 
tate the measurement bases. For the MBQC search al- 
gorithm we require only the discrete set of rotations of 
S = {0,7r}. This implies we are measuring each qubit in 
the basis \± 3 ) = (|0) ± e is \\))/y/2, i.e. the standard |±) 
basis for S = or the |=p) basis for S = tt. 

The MBQC search proceeds by implementing sets of 
measurements (iterators) on m vertical columns each 
consisting of n qubits, on a general nxm cluster state. 
The iterative measurements act sequentially from left to 
right on the (m — 1) columns, leaving the last (i.e. right- 
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most) n x 1 column in the output state. Since the single- 
qubit measurements can be performed simultaneously, we 
consider the set of measurements on a single column of 
qubits to be one operation (iteration) in terms of compu- 
tation time (iteration counts). The number of elements 
2 n in the search list is determined by the "vertical" size 
n of the cluster state grid. The speed of the algorithm 
is completely determined by the "horizontal" length of 
the cluster state grid, or more specifically by the number 
of columns m in the cluster state grid and therefore is 
equal to the the number of sequential iterators. Thus, 
the question of the execution speed of the MBQC search 
is reduced to the identification of the minimum length m 
of the cluster state of width n required to produce a par- 
ticular solution state to the search problem. In this work 
we find that a solution state can be found using m = n 
vertical columns. 

Numerically, each single qubit measurement B(a) is 
a product of a rotation matrix R(a) times a Hadamard 
matrix H, B(a) = R(a) H where 

v ; \isin[a/2) cos{a/2) J v2 \ 1 - 1 / 

(3) 

The unitary matrices that implement the measurement 
based search algorithm takes the following forms for the 
two search sizes shown in Fig.(pQ), 

Ubox = B x {a) (g) B 2 (P) ®H 3 ®H 4 
U C 9 = Bi(a) B 2 (f3) <g> B 3 (7)® (4) 
At (a) ® B 5 (b) <g> B 6 (c) <g> H 7 <g> H 8 <g> H 9 . 

The output state of the calculation is then, 

\*out) = U\* in ), (5) 

a vector of size 2^. \^ ut) can be divided into 2 N ~ n 
sections of size 2 n , each of which correspond to a set 
of specific measurements outcomes on the m columns 
of n qubits. In the first section of 2 n elements of 
l^owt), each of the prior n(m — 1) qubits have been pro- 
jected "correctly" into the state \+ s ) (i.e |+) for s = 0, 
and |— ) for s = tt ) for their given value of s (i.e. 

|*out) = |+si)l ® • • • \+8 N - n )N-n ® \lpout)(N-n)+l,...,N 

where \ipout) (N-n)+i,...,N is the output state on the last 
column of n qubits) yielding the output for the MBQC 
search algorithm. Each of the remaining segments, each 
of size 2 n , correspond to "incorrect results" in which one 
or more of the n(m — 1) qubits has been projected into 
the \— s ) for their given value of s. However, like any 
MBQC algorithm these errors are known and they can 
be corrected during or after the calculation with polyno- 
mial resource overhead [H|, 0] . The essential point is that 
these "incorrect outputs" are calculable and can be triv- 
ially mapped back into the correct output based on the 
measurement outcomes for any size system. Therefore, 
we can assume without loss of generality that the single 



qubits measurements all give the "correct" output. This 
eliminates the need to deal with trivial feed forward cor- 
rections at the end of the algorithm. We can apply the 
projection matrix P = |+ s )(+ s | to each measurement. 
The measurement operations for the 4 qubit box, (j4j is 
then, 




U box = PB 1 (a)®PB 2 (P)®H 3 ®H A , (6) 



This projective operation populates the unitary matrix 
Ubox with the non-zero values in the first 2 n rows and 
zeros in the rest of the matrix. This is convenient for 
calculations involving large cluster states for which the 
use of dense matrices of size 2^ x 2^ is computational 
expensive. We can simply assume the measurements in 
the MBQC scheme proceed "correctly" and truncate the 
matrix to size 2 n x 2^. Note that U is unitary in the 
absence of the projection matrix P. We sample the 2 n 
output by simply rotating the output qubits, the last n 
qubits, to the |d=) basis. This is equivalent to preforming 
a simple Hadamard on each output qubit and measuring. 

Tagging by the Oracle - In GSA the tagging operation, 
which is preformed by the oracle, is a unitary transform 
which flips the phase of a the solution element. Here 
we consider searching for a unique element in a 2 3 = 8 
element list. This tagging operation forms half of the 
Grover iterator. The tagging is one to one, i.e. there 
is a single unique tag for each and every element in the 
search space. In order for the inversion about the mean 
to take place the oracle must continue to apply the same 
tagging operation at each step (iteration) of the search. 

The MBQC search is unique only up to overall phases 
such as ± signs and factors of z, and the final output 
states can be degenerate. For the 8 element search on a 
3x3 grid of qubits, the the final possible measurement 
outcomes on the last column of 3 qubits are 8 fold de- 
generate (ignoring overall phases). In the special case 
of the 4 element search implemented on a 2 x 2 grid of 
qubits, the 4 tagging operations are one-to-one with 4 
search elements. We have found that a one-to-one prop- 
erty associating a tagging operation to a unique output 
state does not scale with the size of the cluster state. 

The MBQC search can be implemented successfully 
for any n with a constant (i.e. fixed) iterator [7]. The 
(n = 3) 2 n = 8 element search on 9 qubits arranged on 
a 3 x 3 grid can be implemented with an iterator com- 
posed of three measurements a,b,c G {0, tt}. This iterator 
would be applied twice to the two columns of 3 qubits, i.e. 
(a,b,c ; a,b,c) on qubits (1,2,3 ; 4,5,6) respectively. How- 
ever, this is not a necessary condition, and the MBQC 
search could instead be preformed with an iterator that 
varies between measurement iterations [7[. 

From Raussendorf, Browne and Briegel [8] we learn 
that the MBQC algorithm can be broken into a series 
of projective measurements and entangling operations. 
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This can be seen in the following relation, 

(P 2 D 2 )(P 1 D 1 ) = (P 2 P 1 )(D 2 D 1 ). 



(7) 



Here {Pi} is a set of n qubit projective measurements, 
thus removing a column of qubit s, and {Di} is a set of 
entangling operations [8] acting on the logical state and 
adding new qubits to the cluster state. The key point to 
note is that the D^s commute with the iVs. 

We used this feature in our calculation of larger clus- 
ter state such as the 81 qubit cluster state (512 element 
search). 
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FIG. 2: Schematic of piecewise cluster state implemen- 
tation. Cluster state map and effective pure state vector in 
MBQC algorithm (see text for description). 

In the (n = 3) 2 3 element MBQC search implemented 
on a 3 x 3 grid of qubits illustrated in Fig. (j2]), the first 
three qubits 1,2,3 are initialized into the |+) state yield- 
ing the unbiased state \ipi n it) Subsequently, the qubits are 
CZ-ent angled with each other, yielding the pure state 
|^o) = CZ lj2 CZ 2j 3\^init). |^o) is the initial system 
pure state of the MBQC search algorithm at iteration 
k = 0, and is the three qubit linear cluster state \C3) Q- 
The MBQC search iterator then begins its first iteration 
(k = 1) by effectively entangling \ipo) with qubits 4,5,6, 
which are also begin in the |+) state, and are themselves 
CZ-entangled with each other, i.e. D\. This creates 
a 2 x n cluster state, similar to the butterfly state |9j. 
The qubits 1,2,3 are then projectively measured, i.e. P\. 
This leaves a three qubit state on 4,5,6 which we de- 
note as l^i), the three qubit system pure state after one 
iteration (k = 1) of the algorithm. In the next itera- 
tion (k = 2) the state is entangled with three new 
qubits 7,8,9, identically to the previous step. For clar- 
ity we call this unitary operation D 2 . It is important 
to note that qubits 4,5,6 are then measured, P 2l and 
the logical pure state \^ 2 } is contained in qubits 7,8,9. 
The projective operations Pi are simply the three qubit 
projective iterators discussed above, and the Di are the 
piecewise construction, i.e. entanglement, of the cluster 
state shown in Fig. ([2]). Qubits 7,8,9 need not exist when 
qubits 1,2,3 are measured since the quantum information 
is entirely stored in 4,5,6. This has particular relevance 
to cluster states constructed from "flying" qubits. The 
significance is that no more than 2n qubits need ever be 
present simultaneously. This constitutes a significant ad- 
vantage for the photon-based MBQC architectures, espe- 
cially since the creation of large photonic cluster states is 
currently limited by the difficulty inherent in producing 
larger numbers of on-demand single photons. 



Algorithm Speed - The speed of the classical search on 
an unsorted list is known to be Q(2 n ~ 1 ). Grover's al- 
gorithm improves upon this lower bound by performing 
amplitude amplification. Grover's algorithm is in essence 
a vector rotation from an unbiased state to a desired tar- 
get state and scales as f£(\/2™) [l|. This quadratic im- 
provement in speed was shown to be optimal for a gener- 
alized amplitude amplification method by Zalka [3]. The 
MBQC search algorithm does not apply amplitude am- 
plification since it utilizes non-reversible projective mea- 
surements Q- 

We now need to determine the "length" of the clus- 
ter state m in terms of the "width" n. To date, there 
exists only one experimental (photonic) implementation 
of the MBQC search algorithm. This was performed by 
Zeilinger et at. [4] for the trivial n = 2, 4 element search 
utilizing the box cluster state in Fig.(pji), i.e. m = n = 2. 
For larger searches n > 3 there is in principle no up- 
per limit to the length of the computation m. However, 
our numerical simulations find that this is not the case. 
Instead, we find that for any width n the cluster state 
must be exactly n columns long. The MBQC search re- 
quires square cluster states; n — 1 measurement itera- 
tions (columns) acting on exactly n qubits at each itera- 
tion, and one (last) column for the output measurement, 
yielding m = n. This result has been numerically tested 
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FIG. 3: Number of iterations (Color Online). The num- 
ber of elements in an unsorted vector and the number of itera- 
tions required for three search algorithms. The classical 2 n_1 
search (straight green line), the standard grover's algorithm 
estimate of L(vr/4)v / 2^J (red line) and the MBQC search algo- 
rithms Log2 2 n (black circle). Searches which are not powers 
of N ^ 2 n are implemented on the next larger cluster state, 
N < 2 n by ignoring the unneeded elements. This is what 
leads to the step like shape of the MBQC curve. 

for 2 n elements from n = 2 —> n = 9. The (n = 9) 
2 9 = 512 element MBQC search requires 81 total qubits, 
72 single qubit measurements and up to 72 post rota- 
tions. Such a large cluster state is currently physically 
unrealizable. Our numerical results show that this clus- 
ter state takes only n — 1 = 8 applications of the MBQC 
iterator. Grover's method requires [(tt/^V^J « 17 ap- 
plications of the amplitude amplification iterator (where 
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\x\ = floor (x)). This result is reflected in figure (j3j) 
where the abscissa has been truncated to maintain clarity 
for the more interesting smaller element sizes. 

Two questions remain; why does the MBQC search not 
give a solution for m < n — 1 iterations, and what hap- 
pens when m > n? To empirically answer the first we 
return to the observation that the MBQC search cannot 
be an amplitude amplification method [H, [?J . Our empir- 
ical results suggest that the MBQC search is more akin 
to a binary decision or search tree. The oracle directs the 
search through the binary tree (from root to leaf) based 
on the oracles application of the tagging operation [7j. 
Each single qubit measurement adds a layer to the bi- 
nary tree. In the n = 2 case there are 2 measured qubit s 
giving 4 possible measurements on qubits 1 and 2 respec- 
tively (0,0); (0,7r); (7r,0); (7r,7r). As both qubits are in 
the same column we can measure them simultaneously. 
The n > 3 case is more complex, one might assume that 
a 3 x 2 (width x length) qubit cluster state, the so called 
butterfly cluster state [9[, might be sufficient. We have 
2 3 = 8 target states and one might suspect that only 3 
"decisions" (single-qubit measurements) are required to 
find the tagged state logarithmically. This is not the case. 
We have tested the 3x2 cluster state numerically and 
find that the output state of such a cluster has uniform 
amplitude for any tag, and is therefore effectively ran- 
dom. This behavior holds for any n and m if the cluster 
state is shorter than it is wide, m < n. 

The 3x3 qubit cluster state (n = 3) from Fig. ([TJd) 
does give a single target output state with probability 1. 
This is somewhat surprising as we have n{n — 1) = 6 "bi- 
nary" measurements to select a solution amongst 2 3 = 8 
distinct solution elements. In this case the 2 6 = 8 x 8 
possible output states are exactly 8 fold degenerate (up 
to overall phase). In addition we find that to obtain a 
solution element, one only needs n of the n(n — 1) mea- 
surements to be variable, i.e. take values S G {0,7r}, 
with the remaining fixed at S = [7]. In this sense, we 
have 2 n distinct tags selecting between 2 n distinct solu- 
tion states, similar to the trivial n=2 case. We find all of 
these discussed behaviors scale with n, at to at least the 
upper limit n = 9 of our computational investigations. 



When m is larger than n we find an interesting cyclic 
behavior. Unlike GSA which cycles between the target 
state and an unbiased state, we find that the MBQC 
search behavior varies significantly. For a given tag, the 
desired MBQC target state will be produced at iteration 
k = n — 1. Afterwards, for some tags, the desired target 
state will be produced every n + 1 iterations. For other 
tags, the states cycle every n + 1 iterations between the 
target state and a second distinct valid output state. Ta- 
ble (|T]) shows the chosen iterator for the case n = 3 and 
the resulting tagged element from 1 to 8, including phases 
at the output of each cycle. Thus each column requires a 
different length cluster state of 3x (number of iterations), 
plus one output column. In Table (pp) the iterators (0,0,0) 
and (7r,0,7r), corresponding to target states 1 and 6 re- 
spectively, give output elements that are always the same. 
The remaining tags cycle between between two n-qubit 
separable target states. The precise cause of this cyclic 
behavior is a topic warranting further investigation. 

Conclusion - We have performed numerical investiga- 
tions of a search on an unsorted list of 2 n elements using 
a MBQC approach that generalizes the experimentally 
confirmed 4 element search performed by Zeilinger et al. 
[4j . Our results lead to the empirical conclusion that the 
MBQC search algorithm is not Grover's Algorithm. More 
significantly, the speed of the two search algorithm are 
very different, with the MBQC search providing an ex- 
ponential speedup over Grover's Algorithm. That is, for 
a 2 n element search, Grover's search algorithm requires 
S7(2 n / 2 ) applications of the Grover iterator acting on n- 
qubits. For the MBQC algorithm acting on n 2 qubits, 
only n—l, i.e. Q(n) applications of the MBQC iterator 
are required, acting on n-qubits at each iteration. 
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are gratefully for useful discussions with M. Fanto, D. 
B. Uskov, and J.R. McDonald. Any opinions, findings 
and conclusions or recommendations expressed in this 
material are those of the author (s) and do not necessarily 
reflect the views of AFRL. 



TABLE I: The tagged search element in the decimal basis and 
the number of iterations at which it is reached. Column I2 is 
the same as Fig. {Dd) and Fig. j2]). 
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